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CN ■ Abstract 

We discuss which part of information about hadron structure encoded in the 
Generahzed Parton Distributions (GPDs) [part of total GPD image] can be restored 
0^ I from the known amphtude of a hard exclusive process. The physics content of this 

partial image is analyzed. Among other things, we show that this partial image 
contains direct information about how the target hadron responses to the (string) 
quark-antiquark operator of arbitrary spin J. Explicit equations relating physics 
^ content of the partial image of GPDs directly to the data are derived. Also some 

Qj • new results concerning the dual parametrization of GPDs are presented. 

(N 
> 

O . Measurements of hard exclusive processes can provide us with rich information about 
hadron structure encoded in generalized parton distributions (GPDs) [1] (for recent re- 
views of GPDs see Refs. [2, 3, 4]). Frequently one speaks about "hadron tomography". 
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Introduction 



O I "3D images of proton", or "fempto- hologram" , etc. Indeed, if we would know GPD as 
^ ■ a functions of all its variable we would determine, among many other things, the dis- 
tributions of hadron's constituents (quarks and gluons) in three dimensional space. The 
amplitudes of hard exclusive processes, which we consider as direct observables, are given 
^ i by the convolution of GPDs with perturbative kernels. It means that measurements of 
I the amplitudes provide us with some kind of sectional images of GPDs - the convolution 
integral "projects out" one of variables in GPDs. This is a typical problem of tomogra- 
phy which is usually solved with help of Radon transformation [5] . Recently the Radon 
transformation was applied in Ref. [6] in order to restore the double distributions from 
known GPDs. 

In the present paper we address the question: What part of the hadron image pro- 
vided by GPDs can be reconstructed if we know the amplitude of hard exclusive pro- 
cess? We give explicit formulae which relate the forward-like parton distributions of dual 
parametrization of GPDs [9] to the amplitude of hard exclusive processes. We shall see 
that at fixed renormalization scale it is impossible to restore the complete GPD image, 
however, part of the total image carries valuable information about hadron structure. 
Furthermore, studying of the scaling violation in hard exclusive processes would allow, 
in principle, to restore total GPD. We suggest a practical way to determine GPDs from 
data on hard exclusive processes. 
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We restrict ourselves to hard exclusive processes at the leading order (LO) of pertur- 
bative expansion as this order provides the bulk information about the hadron structure. 
The perturbative corrections can be included into our analysis. Recently similar problems, 
from other point of view, were addressed in Refs. [7, 8] including the NLO correction. 



Basics of the dual parametrization of GPDs 



For our analysis we employ dual parametrization of GPDs suggested in Ref. [9]. This 
parametrization is based on representation of parton distributions as an infinite series of 
t-channel exchanges [10]. The dual parametrization has been already used to describe the 
data on deeply virtual Compton scattering in Ref. [11]. In this section we remind the 
basics of the dual parametrization and present also some new results. 

In the dual parametrization GPDs H{x,^,t) is expressed in terms of set of functions 
Q2Ax,t) = 0,1,2,...) as [9]: 
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where the function H^'^\x,^,t) is defined on the interval — < a; < 1 and it is given as 
the following integral transformation: 
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The resulting GPD automatically satisfy polynominality property, all sum rules and have 
correct scale dependence. 

Let us remind the main properties of the forward-like functions Q2u{x,t) *: 

• At the LO scale dependence of functions Q2u{x, t) is given by the standard DGLAP 
evolution equation, so that these functions behave as usual parton distributions 
under QCD evolution. 

• The function Qo{x, t) at t — is related to the forward distribution q(x) as: 



Qo{x,t = Q) 



X dz 

{q{x) + q{x)) - - —{q{z) + q{z)) 

^ Jx Z 



The inverse transformation reads: 



q{x) + -q{x) = Qo{x) + ^ Q,{y) . 



(4) 



(5) 



• The Mellin moments of functions Q2v{x^t) give the form factors of the twist-2 local 
operators with fixed conformal spin, see details in Ref. [9]. In particular, the x- 
moment of the lowest functions QQ{x,t) and Q2ix.1t) gives the form factors of the 
quark energy-momentum tensor: 



['dxxQo{x,t) = \M§{t), 
Jo b 

/ dxxQ2(x,t) = di(t) + — M^(t). (6) 
Jo 12 

Here M^(t), di{t) are form factors of the quark part of the energy momentum tensor 
in notations of Ref. [13]. M;^(0) is the fraction of the hadron momentum carried by 
quarks, di{t) is the leading coefficient in the Gegenbauer expansion of the D-term 
[12] and this form factor can be related to the spatial distribution of strong forces 
inside the hadron [13]. In the case of spin-1/2 hadron the above sum rule contains 
one additional form factor J'^(i), with J'5(0)-the angular momentum in the hadron 
carried by the quarks. 

• The coefficients dn{t) of the Gegenbauer expansion of the D-term 

CXD 

D{z,t) = {l-z')Y: d^^)C^/'\z), (7) 

n = l 

odd 

can be computed with help of the following generating function: 

^ d^it) = - / - E(« ^rQ2.iz, t) , ^ , , - 5.0 (8) 

^ a Jo z \^/l + a^z^ J 

odd 

*We restrict ourselves to the case the spinless hadron and singlet distributions. Generalization for 
other cases is trivial 
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Now we give couple of new results concerning the dual parametrization. Using the ex- 
pression (1) we can write the expansion of the GPD H(x,$,,t) around the point ^ = 
with fixed x {x > 0) to the order as follows: 



H{x, e, t) ~ ^ Qo{x, t) + ^ J' ^ Qo{y, t) 



+ 
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It is very instructive expression, we see that the expansion of the GPD in small C, to the 
order ^'^'^ involves only finite number of functions Q2fi{x,t) with fi < u (to the order 
these are only Qo and Q2)- If the expansion of GPD around ^ = is known (for example 
if the GPD is computed in a model), than it can be used to determine the functions 
Q2v{x,t) performing the small ^ expansion order by order. This justifies the naming the 
functions Q2u{x^t)-ioiwaid-\\ke parton distributions. 

Another limiting case, when the Eq. (2) simplifies is the limit ^ —> 1. In this limit 
GPD has all properties of distribution amplitudes. This limit can be used to constrain 
functions Q2u{x,t) because sometimes one can derive the form of the GPD in this limit, 
for example pion GPDs are constrained by the soft pion theorems [10]. 



H{x,i = l,t) = - (l-a;2)E / dy 
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Expanding the integrand in the above equation in variable y we obtain the Gegenbauer 
series for H{x,^ — l,t): 



H{x,^ = l,t) = {l-x')J2h,{t)Cji 



(3/2), 
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with Gegenbauer coefficients computed as: 
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In particular: 
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hi{t) = f'dxx {Qo{x, t) + Q2{x, t)) = ^ M^{t) + di{t) , 
Jo 4 



(13) 



i.e. this coefficient is fixed completely in terms of form factors of the energy-momentum 
tensor. This, actually, is not surprising as it is a consequence of the polynomiality condi- 
tion. 

Finally, we note that to ensure the finiteness of the D-term and existence of the Mellin 
moments of GPDs the small x behaviour of the function Q2u{x,t) ~ should be with 
X<2iy (for i/ > 1) and A < 2 for i/ = 0. 



Forward-like distributions in terms of amplitudes: the inversion problem 



The knowledge of GPD as a function of three variables is equivalent to the knowledge of 
set of forward like parton distributions Q2u{x,t). In this section we attempt to express 
the forward-like parton distributions directly in terms of the amplitudes of hard exclusive 
processes. 

The leading order amplitude of hard exclusive reactions is expressed in terms of the 
following elementary amplitude^: 



A{i,t) = ['dxH{x,C,t) 
Jo 



^-x-iO ^ + x-iO 



(14) 



We see that the amplitude is given by the convolution integral in which dependence of 
GPDs on variable x is "integrated out" . We shall see below, that mathematically from 
the equations (14) one can not completely restore^ the GPD H{x, ^, t). So we are not able 
to perform "complete imaging" of the target hadron from the knowledge of the amplitude. 
The key question is: what part of the "complete image" can be restored from the known 
amplitude? What is the physics content of the restorable part of the complete image? 
Below we address this question in terms of dual parametrization. 

We can express the amplitudes in terms of forward-like functions Q2u{x, t) as following 

[9]: 
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^We restrict ourselves to the singlet (even signature) amplitudes, generalization for odd signature 
amplitudes is trivial. 

■I-We consider the inversion problem at a single value of the large photon virtuality. 
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Here we introduced the function: 

oo 

N{x,t) = Y^x^''Q2u{x,t), (16) 

and the D-form factor: 

D{t)^f:d^{t)^H" dz^. (17) 

n=l ^ ''-^ ^ ^ 

Here D{z,t) is the D-term [12]. 

Now we clearly see that the knowledge of the LO amplitude is equivalent to the 
knowledge of the function N{x,t) and D-form factor D{t). Moreover the D-form factor 
can be computed in terms of N{x, t) and Qo{x, t). Note that the latter function is to great 
extend is fixed by the forward parton distributions, see Eq. (4). Indeed, if we use Eq. (8) 
at q; = 1, we can write: 

D{t) = f - Qo{z,t) (^=L= - lU /' ^ [N{z,t)-Qo{z,t)] -j2== (18) 
Jo z \\/l + z^ J Jo z a/I + ^ 

Note that in this equation all integrals are convergent if the functions Q2u{x, t) satisfy the 
small X behaviour discussed at the end of the second section. In Ref. [16] the following 
representation for the D-form factor in terms of GPD has been suggested: 

/I fix 
— [H{x,x,t)-H{x,0,t)] . (19) 
-1 X 

Unfortunately, this expression is divergent, as H{x,x) and H{x,0) have different coef- 
ficients in front of leading small-x asymptotic, see Refs. [14, 9]. The reason for this 
divergency is that the above expression was obtained by small ^ expansion of GPDs, 
which explodes if x ~ ^. 

Another remarkable physics feature of the function N{x, t) is that its Mellin moments 
are related to the contributions of states with fixed angular momentum in the i-channel: 

dx x'-^ Nix, t) = \ C dz^^^ , (20) 
^0 2 J-i 1 — z 

where $j(^, t) is the distribution amplitude corresponding to two quark exchange in the t- 
channcl with fixed angular momentum J. The quantity on RHS of Eq. (20) carries valuable 
information about the hadron structure - it tells how the target nucleon responses to the 
well defined quark-antiquark probe of arbitrary spin J. 

Now we turn to the inversion problem: how to obtain the function N{x, t) if we know 
the amplitude of a hard exclusive process. This problem is the central for the physics of 
hard exclusive processes. To solve the problem we start with the expression (15) for the 
imaginary part of the amplitude. It is useful to make the following substitution for the 
integration variable x: 
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This substitution correspond to famous Joukowski conformal map [15], which historically 
was used to understand some principles of aerofoil design. After this change of variables 
the expression for the imaginary part of the amplitude gets the form: 



Im^(e,t) = / — MKt)^_, (21) 
^ V ^ 



where the function M(w, t) is related to the function N{x, t) as: 



V ^ / ^2(1-^2)^1- vr^^ 

Now we can easily invert Eq. (21), this can be done in many ways. One of them is to 
view the problem as the tomography problem and solve it with help of inverse Radon 
transformation [5]. The solution of the integral equation (21) is: 

Performing trivial change of variables we arrive at the final result for the function N[x, t): 



This remarkable formula allows to restore the function N{x, t) from the measured imagi- 
nary part of the amplitude. Note that the inversion formula contains the amplitude only 
in the physical region. At ^ — > 1 the imaginary part of the amplitude should go to zero. 
Let us assume that Im A(^,t) ~ (1 — as ^ — > 1, than from Eq. (24) one can easily 
obtain that 

as a; — > 1. 

As to small x behaviour, the Regge like asymptotic of the imaginary part of the 
ampthtude Im A{^,t) ~ 1/C", according to Eq. (24), corresponds to the following small 
X behaviour of N{x, t): 

.r/ N 1 r(l + Q;) 1 



2° T{\)T{\ + a) x^ ' 

Note that the forward-like function Q2u{x,t) enter the function N{x,t) with the weight 
x^'^. Given the same small x behaviour of these functions, the contribution of higher 
functions to the small ^ (large energies) behaviour of the amplitude is suppressed by ^2''. 
This resembles the contribution of daughter Regge trajectories of dual resonance models 
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(DRMs) for hadron interactions. The relations of the dual parametrization of GPDs to 
DRMs will be discussed elsewhere. 

Let us give also the expression for the Mellin moments of the function N{x, t) which, 
as we discussed, are related to the quark exchanges in the t— channel with fixed angular 
momentum, in terms of the amplitude. The corresponding expression has the form: 



I'dxx'-' N{x,t) = ^l^'^ |ilm^(e,t)-e^Im^(e,t)| RAO, (25) 
with the function Rj{C) given by the following integral: 



Jo y/w \ W J V^ — W 

Its small ^ behaviour has the form: 

/eV^^ r(j + i)r(i) 
""'^^^^[2) r(j+|) • 

Finally, we can in principle compute the D-form factor with help of Eq. (18). For 
that calculation we need additional information - the function QQ^x^t). Fortunately this 
function can be fixed (up to t-dependence) with help of Eq. (4). We do not give the 
complete expression for the D-form factor, but rather the key integral in Eq. (18): 



/' - iV(a;, t) = ^ /' f I ilm A{i, t) - ^^Im A{^, t) ] , (28) 

Je X VI + a; Tl J2e i \2 dt, J 

where we introduced infinitesimally small regulating parameter e to cut possible diver- 
gencies at small x. In principle, from the measurements of the imaginary part of the 
amplitude we can determine the function N{x, t) and from measurements of the real part 
the D-form factor D{t). These two measurements would allow to constrain considerably 
the function QQ{x,t) which gives directly the probability density of partons in 3D space 
[19]. 

To summarize this section: all above formulae give us the possibility access the wide 
set of hadrons properties directly in terms of amplitudes of the hard exclusive processes. 
Still we are not able to restore the complete GPD as for this we need to determine not 
only the function N{x,t) but rather the whole set of functions Q2v{x,t). This can be 
achieved if we would know more general function: 



00 

A^a(2:,t) = EM''Q2.(^,i). (29) 

In principle, this function can be determined if we would measure the logarithmic scaling 
violation for the amplitude. Probably, in practice, it is difficult problem and more practical 
way is to model the function Nq,{x, t) using the non-perturbative information about hadron 
structure. We shall discuss possible ways for such modelling after we derive dispersion 
relations and crossing in the next section. 
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Dispersion relations and crossing 



Up to now wc made use of only the imaginary part of the amphtude, let us study the real 
part of the amplitude. For this we take the expression (24) for the N{x,t) in terms of 
lmA{^,t) and substitute it into Eq. (15). After some simple calculations we arrive to the 
following expression for real part of the amplitude: 



Re^(e, t) = 2D{t) + -vp [' dC Imyl(C, t) 

TT Jo 

in which we can immediately recognize the dispersion relation for the amplitude with 
one subtraction at non-physical point ^ = 00 (corresponding to v — {s — u)/Am — 0). 
The D-form factor is the corresponding subtraction constant. This result was obtained 
recently in Refs. [16, 17] by independent methods. We see that the dual parametrization 
automatically ensures the dispersion relations for the amplitudes. 

The very idea of the dual parametrization of GPDs in terms of t-channel exchanges 
was motivated by the crossing relations [10] between GPDs and generalized distribution 
amphtude [18]. The later enter the description of the hard exclusive processes in the 
cross channel, like 7* + 7 ^ /i + In the LO the amplitude of the cross process can be 
expressed in terms of the function N{x,t) analytically continued to time-like t {t > 0)^: 

A-°-(77,i)= l"-N(x,t) 
Jo X 

Here —rj is directly related to the cos(^cm)- cosine of scattering angle in centre of mass 
system, see for details Refs. [18, 10]. Now substituting our inversion formula (24) into 
this expression we obtain, rather simple result: 

A-°-(77, t) = - /'"' d^-^ 1mA (f t] + 2 D{t) . (32) 

TT Jo 1 ~ S V 1^1 / 

Actually this equation is the consequence of the dispersion relation (30). 
A way to model Q2v{x, t) 

As we mentioned above the complete knowledge of GPD is equivalent to the knowledge 
of function (29). In this section we consider possible ways to model this function. For 
simplicity we do not consider the t— dependence and that why we do not write the 
corresponding argument. Without loss of generality we can represent the function as 
the Mellin convolution: 

Q2A^) ^ C- Qo{z) P. (-) , (33) 

Jx z \z J 

§Such continuation can be performed with help of dispersion relations in variable t, see e.g. Ref. [10] 
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then the function Na{x) can be written as the integral transform of the Qo{x): 




where the integral kernel is: 

K{x,y) = ^x'^P,{y). (35) 

1^=0 

The presentation of Q21/ as the Mellin convolution (33) has an advantage as the functions 
Pi,{y) and hence the integral kernel (35) are scale independent at LO. It is little known 
about the kernel (35), its small x expansion has the form: 



K{x,y)r^5{l-y)+x' Pi{y) + ... . (36) 

We make the following ansatz for the integral kernel K[x,y): 

K{x, y) = (1 + A{x))S{l -y) + B{x) . (37) 

Surely, this ansatz is very simple, but still it gives enough freedom in modelling of the 
function Na{x, t). Let us assume that apart from the knowledge of the amplitude we have 
an additional non-perturbative information. For instance, one possess the knowledge of 
the D-term. Using the generating function (8) we easily obtain its relation to functions 
A{x) and B{x) in our ansatz (37): 

^ „ , 1 rfx _ , , f 1 + A{ax) ^ dz B{az) \ 

^ a Jo X [ V 1 + Ci^x^ JO z y/l + a^z^ J 

If the functions A{x) and B{x) are related to each other by the following equation: 



(1 + x''){xA'{x) + B{x)) = x'(l + A{x)) , (39) 

the corresponding D-term is zero. This relation can be used for modelling of the GPDs 
for which the D-term is identically zero, e.g. the combination H{x, {) -|- E{x, {) of the 
nucleon GPDs. For general case we can use the following ansatz: 



A{x) = C{x) + A^°\x), B{x) = B^'^\x), (40) 

where A^^\x) and B^^\x) satisfy Eq. (39). The function C{x) can be determined from 
known coefficients (/„ and known function Qo{x) through the following relation: 



n^i a Jo X \V1 
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If we look for the solution of Eq. (39) in class of polynomials of finite order 2N ^ the 
resulting polynomials A^^\x) and B^'^\x) depend on A*" — 1 free parameters. For example 
the solution in the space of the 4th order polynomials depends on one free parameter a 
and has the form: 

{x) = ax^ + (1 - a)x\ {x) = (1 - 2a)x^ - 3(1 - a)x^ . (42) 

The expression for the function Na{x) in this case has the form: 

N!^\x) = [1 + a {axf + (1 - a) (ax)^] Qo{x) 

+ [(1 - 2a) {axf - 3(1 - a) [axf] / — Qo{z) . (43) 

Jx Z 

Allowing solutions in the class of higher order polynomials we introduce new free param- 
eters which can be fitted to presumably known function N{x), the latter can be obtained 
from the data with help of Eq. (24). 

Conclusions 

In the framework of dual parametrization of GPDs we derive explicit equation (24) which 
expresses particular combination (16) of forward-like functions Q2v{x^t) in terms of the 
imaginary part of the process amplitude. The real part of the amplitude contains one 
additional constant - the D-form factor, which is actually the subtraction constant of the 
corresponding dispersion relation. We show that the D- form factor can be computed 
in terms of the function N{x,t) and the forward-like function Qo{x,t) which is at t = 
is related to forward parton distributions by Eq. (4). It means that the amplitudes of 
hard exclusive processes (at fixed hard scale) provide us only with partial image of total 
GPDs. This partial image is usefully encoded in the function N{x, t) and D-form factor 
D{t) [equivalently in N{x,t) and forward-like function QQ{x,t)]. The function N{x,t) 
contains direct information about the two quark exchange amplitude in t-channel with 
fixed angular momentum. We showed that the J's Mcllin moments of N{x, t) gives the 
response of the target hadron to the (string) quark-antiquark operator with spin J. 
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